Quantum transport through mesoscopic disordered interfaces, junctions, and 

multilayers 



Branislav K. Nikolic* 

Department of Physics and Astronomy, SUNY at Stony Brook, Stony Brook, NY 1 1794-3800 

The study explores perpendicular transport through macroscopically inhomogeneous three- 
dimensional disordered conductors using mesoscopic methods (real-space Green function technique 
in a two- probe measuring geometry). The nanoscale samples (containing ~ 1000 atoms) are modeled 
by a tight-binding Hamiltonian on a simple cubic lattice where disorder is introduced in the on- 
site potential energy. I compute the transport properties of: disordered metallic junctions formed 
by concatenating two homogenous samples with different kinds of microscopic disorder, a single 
strongly disordered interface, and multilayers composed of such interfaces and homogeneous layers 
characterized by different strength of the same type of microscopic disorder. This allows us to: 
contrast resistor model (semiclassical) approach with fully quantum description of dirty mesoscopic 
multilayers; study the transmission properties of dirty interfaces (where Schep-Bauer distribution 
of transmission eigenvalues is confirmed for single interface, as well as for the stack of such inter- 
faces that is thinner than the localization length); and elucidate the effect of coupling to ideal leads 
("measuring apparatus") on the conductance of both bulk conductors and dirty interfaces When 
multilayer contains a ballistic layer in between two interfaces, its disorder-averaged conductance 
oscillates as a function of Fermi energy. I also address some fundamental issues in quantum trans- 
port theory — the relationship between Kubo formula in exact state representation and "mesoscopic 
Kubo formula" (which gives the zero-temperature conductance of a finite-size sample attached to 
two semi-infinite ideal leads) is thoroughly reexamined by comparing their answers for both the 
junctions and homogeneous samples. 
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I. INTRODUCTION 

The experimental discovery&i of a giant magnetore- 
sistance (GMR) phenomenon has revived interest in the 
transport properties of macroscopicahW inhomogeneous 
conductors, such as metallic junctionsu or multilayers.u 
Furthermore, unusual systems for traditional transport 
theory, like single dirty interfaceO which are ubiquitous 
elements in such circuits, have entailed the introduction 
of new concepts to replace usual quantities (e.g., mean 
free path £) used to describe transport in bulk sam- 
ples. Although some of these problems were formulated, 
long ago within the (semiclassical) transport theory,ETEl 
new attacks have employed all (quantum and semiclas- 
sical) transport formalisms developed thus far reveal- 
ing that such problems are by no means resolved Ou In 
particular, the reexamination of various fundamental is- 
sues in the transport theory has been brought about by 
the experimental and theoretical advances in mesoscopic 
physicsJij Thus, the Landauer-Biittiker scattering for- 
malismty has been freauently invoked to study transport 
in both . n o ni -magneticLj and magnetic multilayered con- 
ductors.EjO Obviously, the thorough understanding of 
transport properties due to purely multilayer + disorder 
effects is prerequisite for the analysis of more complicated 
phenomena in inhomogeneous structures. 

Besides providing the means to compute the (quan- 
tum) conductance of finite-size samples, mesoscopic 
methods give additional physical insight by delineating 
transmission properties of the sample (note that one can 



technically treat transport in macroscopic samples with- 
out really using such fully quantum-mechanical transport 
theoriesa) . The finite-size of mesoscopic systems play an 
important role in determining the conductance through 
scattering approach, but no further limitations exist — 
the exactness of results obtained in this guise is heavily 
exploited in throughout the paper. Practical realization 
of this program appears in different incarnations, i.e., dif- 
ferent Landauer-typeE.3 or Kubot£l formulas for a finite- 
size phase-coherent sample attached to semi-infinite ideal 
(i.e., disorder-free) leads. These prescriptions are usually 
made computationally efficient .by. combining them with 
some Green function technique. EZH13 

Here I employ mesoscopic quantum transport methods 
to calculate the conductance of disordered samples which 
are macroscopically inhomogeneous, i.e., composed of dif- 
ferent homogeneous conductors ("layers") joined through 
some interfaces ( "monoatomic layers" ) . In homogeneous 
conductors, whose properties are well-studied throughout 
the history of localization theoryJLj impurities generate 
only microscopic inhomogeneity on the scale ~ Xf (Fermi 
wavelength). Our goal is twofold: (1) Most mesoscopic 
studies -have been focused on the bulk homogeneous con- 
ductoral3 in the weak scattering transport regime. Only 
recently the systems like metallic multilayers Ji3 or single 
dirty interfaceaJ'Ej have been tackled in this guise. By 
employing nonperturbative numerical methods we can 
access strongly disordered junctions, single strongly dis- 
ordered interface (when stacked together to form a bulk 
conductor our interfaces would form an Anderson insu- 
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lator), and multilayers composed of such interfaces and 
bulk diffusive or ballistic mesoscopic conductors. A sys- 
tem is called mesoscopic if its size L is smaller than 
the dephasing length L^, which is a typical distance for 
electron to travel without loosing its phase coherence 
(L<f, < 1 /im in current low-temperature experiments), 
and is therefore determined by decoherence processes 
caused by the coupling to environment either through in- 
elastic scattering (electron-electron and electron-phonon) 
or just be the change of the environment quantum state 
(e.g., spin- flip scattering from a magnetic impurity). The 
term metallic implies that conductance G of a bulk ho- 
mogeneous sample is much larger than the conductance 
quantum Gq — 2e 2 /h. For interfaces one needs a differ- 
ent nomenclature: they are termed "dirty"B if their con- 
ductance G/Gq is much smaller than the number of con- 
ducting channels (N c h ~ k 2 F S in three-dimensions, 
with fci? = 27t/Af, and S is the cross section of a sam- 
ple). Lacking a better language, I denote the multilayers 
studied here as "dirty metallic" , meaning that scattering 
is due to a random potential and their conductance is 
Gm <C GqN^. Nevertheless, layer components are cho- 
sen to be metallic G ^ Gq, and are well-described, in 
the [diffusive limit L ^> £, by semiclassical transport the- 
ory.!^ In terms of material parameters, the resistivity of a 
bulk homogeneous material, from which nanoscopic lay- 
ers are cut out, is few 100 /ificm, which is typical of dirty 
transition metal alloys. Thus, the possible small con- 
ductance of layers is not necessarily caused bv_approach- 
ing the localization-delocalization transitioned upon in- 
creasing disorder, as is usual in the homogeneous bulk 
samples. The conductors are modeled by a tight-binding 
Hamiltonian with on-site potential disorder. This corre- 
sponds to a model of free electrons (understood here as 
Landau quasiparticles with parameters renormalized by 
both band structure effects and Fermi liquid interaction) 
with random point scatterersjj which is used frequently 
in studies of similar systems. While isotropic scattering 
sets semiclassical vertex corrections to zero (which are de- 
termined by ladder diagramsEa generating difference be- 
tween momentum relaxation time and elastic mean free 
time, or scattering-in term in the Boltzmann theory), it 
does not eliminate the higher order quantum interfer- 
ence vertex corrections. These terms are nonlocal op. the 
mesoscopic length scale L^, and therefore invalidateo the 
concept of local position dependent conductivity u{x) as 
usual way of describing the multilayered structures (in 
semiclassical approximation) Since I exploit here exact 
techniques, all quantum locahzatioja-offccts which are not 
necessarily small in dirty systemsJj'Ej are included from 
the onset. All three types of samples are studied for 
electron transport perpendicular to the layers (or inter- 
faces), which is the so-called current perpendicular to 
plane (CPP) geometry.a Once the disorder-averaged re- 
sistance of the multilayer is computed, we can_cpmpare 
it to the resistance given by the resistor model,Efl (i.e., a 
sum of the bulk layers and interface resistances connected 
like classical Ohmic resistors in a series). (2) I investigate 



some fundamental issues in the quantum transport the- 
ory using dirty metallic junctions from above as a testing 
ground, as well as homogeneous disordered samples as a 
reference. That is, I compare the transport properties 
computed from the Kubo formula in exact single parti- 
cle state representation (which was widely usedEEl in the 
"premesoscopic" eraE£l of the Anderson localization the- 
ory) and "mesoscopic Kubo formula" foii-the open finite- 
size system attached to two ideal leads.ll3 In the former 
case the system is closed and the eigenproblem of the 
Hamiltonian is solved exactly by numerical diagonaliza- 
tion, while in the latter case energy levels of a sample 
are broadened by the coupling to the leads, and I use 
real-space Green functions for such open sample plugged 
into the Kubo formulaES (which is then equivalent to the 
Landauer two-probe formulaO) to get its conductance. 
Also, the influence of the leads ("measuring apparatus") 
and the lead-sample coupling on the conductance (which 
is akin to the-problems encountered in quantum measure- 
ment theoryEj) is explicitly quantified. 

The paper is organized as follows. The first Sec. || 
deals with some general remarks on the Kubo linear re- 
sponse formalism and current conservation. This should 
serve as a guidance for proper application of Kubo for- 
mulas on the finite-size systems. In Sec. [II, the Kubo 



formula in exact state representation is compared to the 
Kubo formula for the finite-size sample by evaluating 
them for dirty metallic junction, as well as for homo- 
geneous samples. Section IV presents the results on the 
transmission properties, and conductance derived from 
them, of a single dirty interface as well as thin layers 
composed of few of such interfaces. The calculations are 
completed in Sec. [v| by examining different types of mul- 
tilayers, containing disordered and/or ballistic layers con- 
catenated through interfaces from Sec. IV. Section VI 



provides a summary of the technical, emphasizing phys- 
ical insights gained from them. 



II. KUBO FORMULAS AND CURRENT 
CONSERVATION 

The basic global transport property, for small applied 
voltage, is linear conductance (or equivalent resistance 
R = 1/G). The conductance G is defined by the Ohm's 
law 



I = GV, 



(1) 



relating the current / to the voltage drop V. Since this 
is a plausible relation for a linear transport regime, more 
information is contained in the local form of the Ohm's 
law 



j(r,a;) - J dr'<7(r,r';a;)-E(r>). 



(2) 



In what follows, the focus will be on the DC transport 
uj — > 0. Although it is possible to treat E(r) as an exter- 
nally applied electric field and then include the effects 
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of Coulomb interaction between, electrons as a contri- 
bution to the vertex correction^ the usual approach is 
to use the total (local and inhomogeneous) electric field 
E(r) = Ei oc (r) = — V/i(r)/e, which is the sum of external 
field plus the field due. to the charge redistribution from 
the system response. c3 The electrons are then treated as 
independent quasiparticles. The electrochemical poten- 
tial /z in non-equilibrium situations (like transport) is not 
a well-defined quantity, and cap-serve only e.g., to param- 
eterize the carrier populations The relation (g) defines 
the meaning of the nonlocal conductivity tensor (NLCT) 
g(r,r') as the fundamental microscopic quantity in the 
linear response theory. This quantity gives the current 
response at r due to an electric field at r'. The require- 
ment of current conservation in the DC transport 



V-j(r) =0, 



(3) 



coupled to Eq. (|^) and together with appropriate bound-, 
ary conditions at the infinity or possible interfaces ,E2I 
makes it possible to find E(r). The conductance (dis- 
sipated power by the voltage V squared) is then given 
by 

G=^J drE(r).j(r) 

Q 

= _L| drdr'E(r)-g(r,r')-E(r'), (4) 

Q 

where f2 is the sample volume, and at first sight requires 
the knowledge of local electric fields within the sample. 
However, it was shown in the course of recent reexami- 
nation of transport theory,E!l driven by the problems of 
mesoscopics, that current conservation imposes stringent 
requirements on the form of NLCT 



V • <r(r,r') = <r(r,r') • V' = 0, 



(5) 



where absence of magnetic field is assumed to get this 
special case of a more general theorem£3 This allows us 
to use arbitrary electric field factors in Eq. (||) , including 
E(r) 7^ E(r') [a homogenous field-,i? = V/L is usually 
implied in the textbook literaturecia] . 

When a finite-size sample is attached to two semi- 
infinite ideal leads, the condition (0) is sufficient to show 
thaiO 
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dS 1 -q{r,r')-dS 2 , 



(6) 
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by using the divergence theorem to push_the integration 
from the bulk onto the boundary surfaces going through 
the leads and around the disordered sample (the integra- 
tion over this insulating boundary obviously gives zero 
contribution because no current flows out of it). The sur- 
face integration in the two-probe conductance formula (|^) 
is over surfaces S\ and S 2 separating the leads from the 



disordered sample. The vectors dSi and dS 2 are nor- 
mal to the cross sections of the leads, and are directed 
outwards from the region encompassed by the overall sur- 
face (composed of Si , S 2 and insulating boundaries of the 
sample). It is assumed that voltage in one of the leads is 
zero, e.g., [hi, = and (ir = eV. The meaning of current 
conservation and expression (^) is quite transparent — 
current / at a given cross section depends only on the 
voltages in the leads (in experiment one either fixes this 
voltages by a voltage source, or fixes the current by a 
current source) and not on the precise electric field con- 
figurations. The formula ^ can be generalizedcZl to ar- 
bitrary multi-probe geometry, while the volume-averaged 
conductance (|j) is meaningful only for the two-probe 
measurement. The expression ([|) is valid even in the 
presence of interactions, where many-body effects can 
be introduced using Kubo formalism to get a(r,r') mi- 
croscopically. However, this route is tractable and use- 
ful especially in the case of noninteracting quasiparti- 
cle systems. When a sample is attached to ideaLJpads 
(see below), this established rigorous equivalence!!!! be- 
tween two different linear response formulations — Kubo 
and Landauer-BiittikeJlJ (to work out the proof, Si and 
S 2 should be placed far enough into the leads so that all 
evanescent modes from the sample rhave died out and do 
not contribute to the conductancecl). In the scattering 
approach to transport, pioneered by Landauer through 
heuristic and subtle arguments, the conductance of a 
noninteracting system is then expressed in terms of the 
transmissions probabilities between different quantized 
transverse propagating modes defined by the leads as 
asymptotic scattering states [see Eq. (28)]. 

Although Boltzmann formalism can provide semiclassi- 
cal expression for NLCT (which is nonlocal on the scale 
of the sample size because of the classical requirement 
of current conservations), the standard quantum route 
to it is the Kubo linear response theory (KLRT). Lin- 
ear response theory gives a full quantum description of 
transport (i.e., it includes-<quantum interference effects 
in the motion of electronsEj) in non-equilibrium systems 
that are still close enough to equilibrium for vanishingly 
small applied external electric field SB Therefore, the for- 
mulas above involve the equilibrium expectation values 
of the corresponding quantum-mechanical operators, in 
accordance with fluctuation-dissipation theorem under- 
lying KLRT [e.g., it is assumed that j(r) is the expecta- 
tion value of the current density operator in the quantum 
formalism]. While the current is response to the total 
electric field (external + induced), the NLCT is obtained 
as a response to the external field onlyt£lE3 because the 
current induced by the external field is already linear 
in the field. Therefore, in the linear response one does 
not need the corrections due to induced non-equilibrium 
charges (electron-electron interaction is taken in the equi- 
librium, e.g. through renormalized parameters of the 
Landau quasiparticles and the self-consistent screening 
effect on the impurity scattering). 

Since Kubo NLCT is not experimentally measurable 
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quantity, macroscopic conductivity a is obtained by vol- 
ume averaging NLCT through Eq. (Q) where a = G L/S 
for a cubic sample of length L and cross section S (limit 
£1 = LS — > oo is assumed, while keeping the impu- 
rity concentration finite). This usually implies a ho- 
mogeneous E(r) factors, which is justified by the fact 
that Kubo expression for NLCT is divergenceless (as can 
be easily provenEH from its current-current correlation 
formj^a) i.e., it satisfies Eq. (||). Thus, the volume- 
averaged conductivity <r relates the spatially averaged 
current j = J dr j(r)/fi to the spatially- averaged electric 
field, 



j = aE. 



(7) 



For a noninteracting system of fermions this leads to a 
Kubo formula in exact single-particle state representa- 
tion (KFESR) 
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a, a' 



\(a\v x \a')\< 



x5{E a -E F )6{E a ,-E F ), 



(8) 



where the factor of two for spin degeneracy will be ex- 
plicit in all formulas. The delta functions in Eq. (|J) em- 
phasize that conductivity is a Fermi surface (E F is the 
Fermi energy) property at low temperatures [T — ► => 
-df(E)/dE ~ S(E F ~ E), f(E) being the Fermi-Dirac 
distribution function]. Here v x is the x-component of 
the velocity operator and \a) are eigenstates of a single- 
particle Hamiltonian, H\a) — E a \a). The velocity oper- 
ator v is determined through the equation of motion for 
the position operator f 



ihv = ih- 



dr 
i— 

dt 



[*,#], 



(9) 



where H is the Hamiltonian before the application of the 
external electric field (in the spirit of FDT). 

In the general case, conductivity is a tensor, but since 
symmetries are restored after disorder-averaging, one can 
use a — (a xx + a yy + a zz )/3 as the scalar conductiv- 
ity. This is valid only in the case of homogenously dis- 
ordered sample. For example, in our metal junction or 
multilayers <r xx is different from <j yy and <j zz . To get 
the conductivity (as an intensive quantity) in KLRT, the 
thermodynamic limit f2 is always understood (no station- 
ary regime can be reached in a system which in neither 
infinite nor coupled to some thermostat). In the disor- 
der electron physics this also bypasses the ambiguity of 
conductivity which scalesEa with the length of the sys- 
tem (although scaling is unimportantca in the metallic 
regime G ^> Gq). Nonetheless, the computation of trans- 
port properties from exact single particle eigenstates, ob- 
tained by the numerical diagonalization of Hamiltonian 
of a finite-size system, has been frequently employed in 
the development of disordered electron physics EB In fact, 
it is still the standard method for computing the con- 
ductivity of many-body systems by diagonalizing small 
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FIG. 1. A two-dimensional version of the actual 3D model 
used here for a two-probe measuring geometry. Each site 
hosts a single s-orbital which hops to six (or fewer for surface 
atoms) nearest neighbors. The hopping matrix element is t 
(within the sample), £l (within the leads), and tc (coupling 
of the sample to the leads). The leads are semi-infinite and 
connected smoothly at ±oo to macroscopic reservoirs biased 
by the chemical potential difference fj,L — = eV. 

clusters of some lattice fermion Hamiltonian.@ However, 
direct application of the formula (||) leads to a trouble 
since eigenvalues are discrete when the sample is finite 
and isolated. Strictly speaking, the conductivity is then 
a sum of delta functions, and to obtain finite conductivity 
they have to be broadened into functions having a finite 
width larger than the level spacing. Thus, there are two 
numerical tricks which can be used to "circumvent" this 
problem: (i) one can start from the Kubo formula for the 
frequency dependent conductivity 
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„ f{E a )-f(E a ,) stJp 
x z o(E a 



E a , - nw), (io) 



average the result over finite u> values, and finally extrap- 
olatec3 to the static limit w — > 0; (ii) The delta functions 
in Eq. (pf) can be broadened into a Lorentzian 
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5{x) -> 5(x) = - 



(W2) 2 



nx 2 + {7 1 /2) 2 ' 



(11) 



where r\ is the full width at half maximum of the 
Lorentzian. I find that both methods produce similar 
results. The calculations presented below use the broad- 
ened delta function S(x) in the formula for static Kubo 
conductivity. 

The Green operator of a noninteracting Hamiltonian 
(with specified boundary conditions), 
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[E — H± irj}~ 



(12) 
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contains the same information encoded in the single- 
particle wave function. Therefore, using the expansion 
of single-particle Green operator G r ' a — J2 a \ a )( a \/(E — 
E a ± irj) in terms of exact eigenstates, the Kubo formula 
for the macroscopic (volume-averaged) longitudinal cop=. 
ductivity can be recast into the following expressionO 



Ti- 



ft 

6{E -H) = --ImG = ^(G r - G a ) 



v x S(E — H) v x S(E — H) , (13a) 

(13b) 
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When employing this formula one faces-the same prob- 
lem as in the KFESR — in order to defineeHl retarded (r) or 
advanced (a) Green functiop-a small parameter r\ — > + 
requires numerical handlings analogous to the one ( JTl] ) 
involved in the KFESR (||). However, it is possible to de- 
rive another, "mesoscopic Kubo formula" , for the finite-, 
size sample attached to two semi-infinite clean leads.cll 
This formula represent a fully quantum-mechanical ex- 
pression for the zero-temperature conductance 



4e 2 1 / ~ 
G X x = —, — rrr Tr ( hv x lm G hv x Im G 
h Li \ 



(14) 



which, being measurable, is the only meaningful quan- 
tity to discuss in mesoscopics (the alternative is to intro- 
duce conductivity as the NLCT). Namely, in quantum- 
coherent samples (L > L^) local description of transport 
in terms of conductivity breaks down because of the non- 
local correction (e.g., weak localizations) induced on the 
scale of ^> £, which is much larger than the elastic 
mean free path. Thus, mesoscopic samples have to be 
treated as a single coherent unit ("giant molecule"), so 
that its conductance cannot be interpreted as a combina- 
tion G = aA/L of the conductances of its parts (i.e., such 
description becomes applicable only at high enough tem- 
peratures). Although the formula (|l4|) apparently looks 
like Eq. (|l3|), the crucial difference is that the following 
Green operator is plugged there instead 



[E- H - E r ' a y 



(15) 



This change is not innocuous: G r (G a = [G r Y , S a = 
[E r ]t) has acquired a self-energy term from the "inter- 
action" with the leads !El This looks like the self-energy 
term in the many-body Green functions, but this one 
is exactly calculable, as shown below. Thus, the leads 
enforce new boundary conditions. Since Y, r (E) contains 
an imaginary part when energy E belongs to the band 
of a lead, no small parameter r\ is needed to define the 
retarded or advanced Green operator. 

The lattice model of a two-probe measuring geome- 
try, to be used in the subsequent sections for evaluation 
of Eq. (|l4|), is shown in Fig. |l|. The three-dimensional 
(3D) nanocrystal ( "sample" ) is placed between two ideal 
(disorder-free) semi-infinite "leads" which are connected 
smoothly to macroscopic reservoirs at infinity. The elec- 
trochemical potential difference eV = Hl — which 



drives the current,@ is measured between the reservoirs. 
It is assumed that reservoirs inject thermalized electrons 
at electrochemical potential /i£ (from the left) or [lr 
(from the right) into the system. All inelastic scat- 
tering occur in the reservoirs, which therefore ensure 
steacU-.state transport in the central region. Semi-infinite 
leadso are a convenient method to take into account 
electrons entering of leaving the phase-coherent sample 
(i.e., effective Hamiltonian of an open system, H + S r 
in Eq. (fl5|), is non-Hermitian) . This makes it possible 
to bypass explicit modeling of thermodynamics of per- 
fect (i.e., unaffected by the flow of current) macroscopic 
reservoirs which are introduced heuristically in the scat- 
tering approach (for a different interpretation of the role 
of perfectly conducting leads in the derivation of Kubo 
formula for a finite-size system see Ref. ^?j). The leads 
have the same cross section as the sample, which elimi- 
nates scattering induced by the wide to narrow geometry 
at the sample-lead interfaceo The whole system is de- 
scribed by the following Hamiltonian with nearest neigh- 
bor hopping integrals t mn 



H = y^g m |m)(m| 



<m,n> 



tmn|m)(n| 



(16) 



where (r|m) is the s-orbital ifj(r — m) located on site m. 
The site representation, defined by the basis states |m), 
can be interpreted either as a tight-binding description 
of electronic states or as discretization of corresponding 
single-particle Schrodinger equation. The sample is the 
central section with N x N y x N z sites. The hopping in 
the sample set the unit of energy i mn = t. The disorder 
will be introduced by taking e m to be a random vari- 
able. The leads are clean (e m = 0) but, in general, have 
different hopping integrals £ mn = t-^ =/= t. Finally, the 
hopping which couples sample to the leads is t mn — tc- 
Hard wall boundary conditions are set in the y and z di- 
rections. The different hopping integrals introduced here 
are necessary when studying disordered samples to get 
the conductance at Fermi energies throughout the whole 
band extended (compared to the clean case) by disorder, 
i.e., in such calculations one has to set > t. 

In site representation, Green operator G r,a becomes a 
Green function matrix G r ' a (n, m) = (n|G r,a |m). The 
self-energy E r ~ Y7 L + S^j "measures" the coupling of 
the sample to the leads. This causes the broadening of 
initially discrete levels and, thus, the finite lifetime of 
electron in the sample. Since electron which leaves the 
sample does not return phase-coherently, the dephasing 
length is by definition — L x . Each lead generates its 
self-energy term which has nonzero matrix elements only 
on the edge layers of the sample which are adjacent to 
the leads. They are defined as 



m ) = *c5£,_R(ns, m s ), 



(17) 



with g r L fl (ns, mg) being the surface Green function of a 
bare semi-infinite lead between the sites rig and mj in 
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FIG. 2. Diffusivity (D x (Ef)) of a disordered binary alloy 
modeled by a tight-binding Hamiltonian with quenched dis- 
order (ea = —2 and eb = 2) on a lattice 18 x 8 x 10: (a) 
computed using the Kubo formula (^) in terms of the Green 
function for the sample with attached leads; (b) computed 
from the Kubo formula in exact single-particle eigenstate rep- 
resentation (Kj), where the width of the Lorentzian broadened 
delta function is ri — 25A(Ef). Also plotted (c) is the IPR, 
defined in Eq. (E5I), which gives insight (i.e., degree of local- 
ization) into the structure of eigenstates used to evaluate the 
Kubo formula Eq. (^3|). The disorder averaging is performed 
over 50 samples. 

the end atomic layer of the leadEl (adjacent to the corre- 
sponding sites n and m inside the conductor). I provide 
here explicit expression for these self-energy terms in the 
most general case when hopping integrals t ^ tc ¥= £l 
are different in different parts of the setup in Fig. n] 

2 2 

S L,fl( n . m ) = Ni + lN +l sm(k y n y a) sin{k z n z a) 

xT, r (ky, k z ) sm.(k y m y a) sin(k z m z a). (18) 

This expression is obtained by expanding the surface 
Green function 9£a(ns,ms) in terms of exact eigen- 
states of the Hamiltonian of a semi-infinite lead, which 
satisfy particular hard wall transverse boundary condi- 
tions. Here (n, m) is the pair of sites on the surfaces 
inside the sample which are adjacent to the leads L or R, 
in accordance with (|T7j). Since electrons are confined in 
the y and z direction, electronic states in the lead have 
transverse wave function which are labeled by the dis- 
crete quantum numbers. They define subbands and cor- 
responding conducting "channels" [i.e., transverse propa- 
gating modes as the product of transverse wave function 
and Bloch wave in the x-direction, (m\k y ,k z ) (g) \k x )], 
which represent a basis of scattering states for the scat- 
tering matrix of the disordered region, as envisaged in the 
Landauer picture of quantum transport. The self-energy 
E r (fcj,, k z ) of the channel (k y , k z ) is given by 

Z r (k y , k z ) = Jjr ~ iyJitl-Elj , (19) 



FIG. 3. The result of the same computation of (D x (Ef)) 
as in Fig. ^, but for a different type of disorder, diagonal one 
of the strength W = 10, introduced in the sample (the same 
labels apply to both figures). 

for \E%\ < 2<l- I use the shorthand notation E^ = E — 
e(k y ,k z ), where e{k y ,k z ) = 2t\\cos(k y a) + cos(k z a)] is 
the energy of quantized transverse levels in the lead. In 
the opposite case \E%\ > 2iL we get 

S r (fc„ k z ) = (ee - sgnEzy/E* -4^ . (20) 

When the level spacing of the subbands is much smaller 
than Ef, the number of channels at Ef is large and the 
finite-size model can describe a metal. The possibility 
to evaluate exactly these self-energies allows us to avoid 
the inversion of an infinite matrix which would formally 
give the Green function for the whole systemEJ in Fig. [I]. 
The real-space Green function technique described here, 
pioneered by Caroli et aZJi3 long before official inception 
of mesoscopic physics, treats an infinite system with a 
continuous spectrum by evaluating only the Green func- 
tion between the states inside the sample. This offers 
an alternative to handling a finite-size sample through 
periodic boundary conditions, whose discrete spectrum 
prevents a direct evaluation of the Kubo formulas (||) 
and (|l3|). In general, the technical advantage of Green 
function techniques is that they can be used even when 
well-defined asymptotic conducting channels in the scat- 
tering approach to quantum transport are hard to de- 
fined (i©ig., this is the case when leads have complicated 
shape,t3 making it hard to get explicitly the exact asymp- 
totic eigenstates and their eigenspectrum). 

The mesoscopic Kubo formula can be evaluated at any 
(continuous) Fermi energy. At first sight, it seems like 
that these leads to a much greater computational com- 
plexity than in the case of KFESR because of the need 
to find inverse matrix in Eq. (|l^) at each chosen en- 
ergy, and then perform the trace over the site states 
Tr(. . .) = E m ( m l(- ■ -)l m >- Namely, in the KFESR one 
needs to diagonalize Hamiltonian only once and then use 
obtained eigestates to compute matrix elements of ve- 
locity operator in the eigenbasis. However, this is only 
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an apparent difficulty once the conservation of current is 
invoked. Since all Kubo formulas for conductance stem 
from Eq. (jj) , one has to assume (or choose) some electric 
field factors therein. Current conservation, i.e., the fact 
that I has to be the same on each cross section, means 
that conductance is independent of this choice. The min- 
imal space for tracing, when system is described by tight- 
binding Hamiltonian (TBH) in Eq. (|ilf), is obtained by 
taking both factors E(r) and E(r') to be nonzero only 
on two adjacent planes of the lattice. Namely, the ex- 
pectation value of the velocity operator (^|) in the site 
representation is nonzero only between the states resid- 
ing on adjacent planes, 



(ml^ln) = -t mn [m x - n x ) , 



(21) 



for a TBH (JlJ) with nearest-neighbor hopping. The fi- 
nal result is then to be divided by a 2 [instead of L 2 in 
Eq. (0)]. 

Thus, only a block matrix 2N z N y x 2N z N y of the Green 
function for the whole sample has to be computed explic- 
itly (e.g., those elements which connect first two planes 
of the sample in Fig. |l|). Because [E — H — S r ] is a 
band-diagonal matrix of a bandwidth 2N y N z + 1, one 
can shorten the time needed to compute this block of the 
Green function (15) by finding the LU decomposition of 



a band-diagonal matrix, followed by forward-backwacd 
substitution for each needed Green functioB-element.EHl 
This requires similar number of operational as pin the 
more familiar recursive Green function techniqueo 



III. TRANSPORT THROUGH DIRTY METAL 
JUNCTIONS 

This section studies the static DC transport properties 
of a metal junction composed of two disordered conduc- 
tors with different kind of disorder introduced on each 
side of the contact interface. Both conductors are mod- 
eled as a disordered binary alloy (i.e., composed of two 
types of atoms, A and B) using TBH of Eq. @. The 
quenched disorder is simulated by taking the random on- 
site potential such that e m is either ea or eb with equal 
probability. Specifically, I take the lattice 18 x 8 x 10 
on each side of the junction and for the binary disorder: 
Ea = —4, eb = on the left; and ea — 4 and eb = 
on the right. Thus, the junction has naturally a rough 
interfaces because of the random positions of three dif- 
ferent types of atoms around the plane in the middle of 
the junction. 

We commence with a description of the junction, as 
well as homogeneous samples used as a reference, in terms 
of traditional Kubo theory based on the KFESR (||). Al- 
though Ref. [hi] stated as one of the motivations for un- 
dertaking the rigorous derivation of Landauer two-probe 
formula from KLRT, that traditional use of Kubo for- 
mulao was numerically demanding, today's computers 
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FIG. 4. The density of states of: (a) clean metal (W — 0); 
(b) dirty metal (W — 6 on a lattice 15 x 15 x 15 
averaged over 50 disorder configurations), obtained as 
N(E) = 2/fl ~}2 a S(E — E a ) from the exact eigenspectrum 
E a of a closed sample Hamiltonian; (c) dirty metal [W = 6 
on a lattice 10 x 10 x 10 averaged over 50 disorder config- 
urations], obtained from the imaginary part Eq. (^) of the 
Green function ([L5|) of an open sample; (d) is computed in 
the same way as (c) except that smaller lattice 4 x 4 x 4 is 
used. 



are much more powerful, and it is of interest to compare 
this method to a moderr££l mesoscopic use of the Kubo 
formula. Kubo theory permits the discussion of diffusiv- 
ity D a of an eigenstate |a). This quantity is extracted 
directly from the KFESR for Fermi gas of noninteracting 
quasiparticles 



e 



E 



df 
8E a 



D a = e 2 N(E F )D. (22) 



Thus, the quantum diffusivity is given explicitly by the 
following expression 



D x a = nhJ2\(a\v x \a')\ 2 5(E a -E al ). 



(23) 



In the semiclassical regime, D a — > = Vk£f./3. How- 
ever, D a can be used to characterize transport even in 
the nonperturbative regime where semiclassical concepts, 
like mean free path, loose their meaning.Eil 

To simplify calculations, I compute the quantum dif- 
fusivity (D x (Ep)) averaged over the disorder ((...) de- 
notes averaging over an impurity ensemble), as well as 
over a small energy interval. This is an additional trans- 
port information, obviously related to conductivity (|2^), 
which is not usually seen in the literature on disordered 
electron physics.E3 It has been exploited in the trans- 
port physics of glasses (e.g., to-study the thermal con- 
ductivity in amorphous silicone!]). The width r\ of the 
Lorentzian S(E a — E a i) in (|2^) is chosen as some mul- 
tiple of the local average level spacing A(E a ) in a small 
energy interval around the eigenenergy E a . The method 
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FIG. 5. Diffusivity (D x (Ef)) for a metal junction com- 
posed of two disordered binary alloys (modeled by TBH on a 
lattice 36 x 8 x 10): (a) computed using the Kubo formula ( p^| ) 
in terms of the Green function for the sample with attached 
leads; computed from the Kubo formula in exact single par- 
ticle eigenstate representation (^), where the width of the 
Lorentzian broadened delta function is (b) rj = 25A(£'f), (c) 
r\ — 10A(_Bf), and (d) r\ — 5A(Ef). Disorder averaging is 
performed over 50 different samples. 



of computing (D x (Ep)) is as follows: a set of eigenstates 
(the number of eigenstates is equal to the number of lat- 
tice sites N s = N x N y x N z ) is obtained by numerical 
diagonalization; t— > for each eigenstate I compute D^, 
where summation is going over all states \a') "picked" 
by the Lorentzian S(E a — E' a ) (centered on E a ) in an 
energy interval of 3?y around E a ; i— ► finally, I average 
the diffusivities over the disorder and energy interval de- 
fined by the bin of the size AE = 0.0225. The efficient 
way of calculating quantum-mechanical average values of 
some operator, like (a\v x \a') appearing in the definition 
of eigenstate diffusivity (|23|), is to multiply three matri- 



ces a' ■ v x ■ a, where a is a matrix containing eigenvectors 
|a) as columns, and then take modulus squared of each 
matrix element in such product. The number of opera- 
tions in the naive calculation of the expectation values 
of operator, where each of them is calculated separately, 
scales as ~ N^, while in the method presented above it 
scales as ~ Nf (N s x N s is the dimension of operator 
matrix). This procedure becomes a natural choice once 
we understand that it actually transforms the matrix of 
the operator v x from the defining (site) representation 
into the representation of eigenstates \a). The end result 
of the calculation, the average diffusivity (D x ), is related 
to the conductivity through Einstein relation 



e 2 N(E F )(D x (E F )). 



(24) 



Here N(Ep) is the density of states (for both spin com- 
ponents) evaluated at the Fermi level Ep. 

I first calculate (D. X (E F )) for the homogeneously dis- 
ordered sample (with binary disorder ea = — 2, £b = 2) 
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FIG. 6. Local density of states (LDOS) integrated over the 
y and z coordinates inside the metal junction studied in Fig. ^. 
This "planar LDOS" p(m x ,E) is computed from the exact 
eigenstates of the TBH using Eq. (p7|). The result is plotted 
after p{m x ,E) is averaged over several planes along the x-axis 
(the planes used in this procedure are given on each panel). 



modeled on a lattice 18 x 8 x 10, as shown in Fig. g. To get 
an insight into the microscopic features of the eigenstates 
used to evaluate the KFESR (a few eigenstates around 
each Ep determines the transport properties at Ep), this 
figure also plots the inverse participation ratio (IPR)l3 



(25) 



averaged over disorder and energy [\I/ Q (m) = (m|a)]. 
This is the simplest single-number measure of the degree 
of localization (i.e., the bigger the IPR the more localized 
is the state is, e.g. IPR= N s corresponds to a completely 
localized states on one lattice site). The LER can also 
be related to the average return probabilityllj that parti- 
cle, initially launched in a state |m) localized on a lattice 
site m, will return to the same site after a very long time 
(which is determined by its diffusive properties encoded 
in D a ). The second calculation, shown in Fig. |[ is for the 
homogeneous sample described by the standard Ander- 
son model where e m £ [— W/2, W/2] is a uniform random 
variable. This results are to be compared to the reference 
calculation based on the mesoscopic Kubo formula. 

Strictly speaking, the concept of eigenstate does not 
exist in open systems. Such systems can be characterized 
by non-Hermitian Hamiltonians,Ej which do not conserve 
total probability (cf. Sec. [nj). This is a consequence of 
a simple physical fact that an electron stays finite time 
within the sample before escaping into the surrounding 
leads (i.e., initial discrete energy levels are broadened by 
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FIG. 7. Two-probe conductance of a disordered conductor 
(averaged over 50 samples) modeled by the Anderson model 
with diagonal disorder W = 6 on a simple cubic lattice 10 3 . 
The computation is done using the mesoscopic Kubo for- 
mula ([w]) for the finite sample attached to two semi-infinite 
leads characterized by two different values of hopping integral 
<l- Note that the conductance vanishes at \E\ = 6i (band 
edge in a clean sample) when £l = t. 



the coupling to the leads). Nonetheless, we can still getE3 
the density of states (DOS) from the imaginary part of 
the Green function jlq ) 



N(E F ) = V ImG r (m,m;B f ). 



(26) 



That such DOS of an open 3D system is indistinguishable 
from the one computed from the distribution of energy 
eigenvalues, N(E F ) = (2/fi) (J2 a S(E F — E a )), is shown 
in Fig. ||. This is in spite of the fact that leads strongly 
perturb the edges of the system, which is clearly exhib- 
ited only on the smallest lattice in Fig. |]. Thus, the aver- 
age diffusivity (D x (Ep)) can be computed in a straight- 
forward manner from the Einstein relation ( p4| ) where 
conductivity is formally expressed from the disorder- 
averaged conductance a xx — (G XX }L X / S given by the 
mesoscopic Kubo formula. 

In both calculations for the homogeneous samples it 
appears that discrepancy between the Kubo formula in 
single particle representation (|^) and the exact method, 
based on the formula ( |l4| ) for the sample + lead system, 
is only numerical. In fact, the difference is very small 
in the disordered binary alloy and a bit larger in the 
Anderson model with continuous disorder. It originates 
from the ambiguity in using the width rj of the broadened 
delta function. Namely, nonzero 77 effectively introduces, 
inelastic scattering as an uncorrelated random event £3 
The increase of the diffusivity close to the band edges of 
diagonally disordered Anderson model (cf. Fig. |^) was 
seen long time ago in the direct simulations of the wave 
functiop-|diffusion performed in early days of localization 
theory.Ea 
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FIG. 8. Two-probe conductance (averaged over 50 sam- 
ples) of the same dirty metallic junction sample as in Fig. ^, 
but attached to two semi-infinite ideal leads with several dif- 
ferent hopping integral tt or lead-sample couplings tc (hop- 
ping integral in the sample sets the unit of energy t = 1). 
The computation is based on the mesoscopic Kubo formula 
Eq. (113). 



The same analysis is repeated for a junction (intro- 
duced at the beginning of this section) which is com- 
posed of two different disordered binary alloys on each 
side. The result is shown in Fig. ^. Large fluctuations of 
the diffusivity (i.e., quantum conductance G from which 
diffusivity is computed at specific Ep) are caused by the 
conductance being of the order of 2e 2 /h (Fig. |J). Such 
fluctuations are less obvious in the result based on the 
KFESR because of the extra averaging over energy pro- 
vided by the Lorentzian broadened delta function. Here 
the discrepancy between the two different methods is not 
only quantitative, but the KFESR (||) shows nonzero 
diffusivity (and thereby conductivity, since global DOS 
is nonzero all throughout the band) at Fermi energies 
at which there are no states on one side of the junc- 
tion which can carry the current .L3 The result persist 
even when the width 77 of the Lorentzian broadened delta 
function is decreased. Therefore, it is not an artifact of 
the numerical trick used to evaluate KFESR. The states 
which have nonzero amplitude all throughout the junc- 
tion cease to exist at \Eb\ ~ 4.7, which is seen by inspect- 
ing the local density of states (LDOS), integrated over y 
and z coordinates 



p{m x ,E) = P( m > E ) = 

m y ,m z 



Such "planar LDOS" along the a;-axis is plotted in Fig. ||. 
It changes abruptly while going from one side of the junc- 
tion to the other side (except for the small tails near the 
interface). On the other hand, Fig. [5] explicitly demon- 
strates that Kubo formula ( |l4| ) for an open finite-size 







sample plugged between ideal semi-infinite leads cor- 
rectly describes this junction. Namely, the diffusivity 
vanishes at the same point at which LDOS goes to zero. 
It should be emphasized that once the leads are attached, 
two new interfaces (lead-sample) in the problem arise. 
The Landauer-Biittiker scattering approach to transport 
intrinsically takes care of these boundaries by consider- 
ing the realistic finite-size system where electrons can 
leave or enter through some surfacesJnto the rest of the 
circuit represented by the leadso£3 Thus, mesoscopic 
developments have clarified the way of proper applica- 
tion of Kubo formulism to finite-size samples (which is 
ultimately related to the eternal puzzle of the origins of 
dissipation in conservative systems, technically the only 
ones appearing in the Kubo analysis). E3 This in turn has 
justified heuristic arguments of Landauer on a rigorous 
basis.BEll 

The conductance of an infinite (sample + leads) system 
will go to zero at the band edge of a clean lead \E^\ — 6t if 
we use the same hopping parameter in the lead t^ = t as 
in the disordered sample. This stems from the fact that 
are no states in the leads which can carry the current for 
Fermi energies \Ep\ > 6t outside of the clean TBH band. 
Technically, the self energy E r is real at these energies, 
which leads to Im G in ( |i"4| ) being zero. Thus, the conduc- 
tance of the whole band of disordered sample cannot be 
computed unless we increase t^ in the leads. This is illus- 
trated in Fig. ^ for the homogeneous sample described by 
the Anderson model where band edge of the disordered 
sample \Eb\ > 6t lies outside of the clean metal band 
defined by crystal symmetry. Thus, a natural question 
arises when using £l, 7^ t'- how sensitive is the con- 
ductance to the properties of leads or the sample-lead 
coupling? This problem resembles the quantum mea- 
surement problem because semi-infinite leads can also be 
viewed as-ar-aiacroscopic apparatus necessary for mea- 
suremenjJ23'E3 This is further stressed in multi-probe ge- 
ometries!!^ where extra leads are introduced to measure 
the voltages along the sample (besides the two leads, used 
here, through which the current is fed). The measured 
conductance is then property of both the sample and the 
measurement geometry, which is one of the reasons for 
mesoscopic physics to employ only sample- and measure- 
ment geometry-specific quantities, like quantum condrusa 
tance, instead of intensive quantity like conductivity.^ 
From a transport point of view, it is clear that lead- 
sample interfaces introduce extra scattering for different 
effective masses or Fermi velocities in the sample and the 
lead. Thus, the exactness of the conductance calculated 
here is, in fact, a feature of the whole sample + leads 
setup in Fig. [l] (that is akin to any kind of measurement 
in quantum mechanics), and it is important to confirm 
that general conclusions about the relationship between 
different Kubo formulas do not depend on particular val- 
ues of the chose parameters for the leads. 

It is understoodEll that if broadening of the energy lev- 
els due to the coupling to semi-infinite leads is greater 



than the Thouless energy E Th ~ HV/L 2 , (V = v F £/3 
is the classical diffusion constant), then level discrete- 
ness of an isolated sample is unimportant and G will 
be independent of the properties of leads. This limit 
corresponds to an "intrinsic" (Thouless) conductance 
G/Gq — 27ri?Tii/A (A is the mean level spacing) that 
is smaller than the conductance determined by the lead- 
samole contact O It requires disordered enough sam- 
pled and. leads of the same cross section as that of the 
sample.Ej Thus, even though attaching the sample to 
the leads brings dramatic changes (discrete spectrum is 
changed into continuous and new boundary conditions 
are introduced), the conductance is determined by the 
sample properties only, and dissipation in the reservoirs 
does not enter into computational algorithm. This de- 
pendence is studied in Fig. ^| by looking at the conduc- 
tance of our model junction as a function of the hopping 
in the leads £l and the coupling tc- The conductance is 
virtually independent of , which is a consequence of the 
smallncss of the disordered junction conductance (as dis- 
cussed above). It goes down drastically with decreasing 
of the coupling tc (the same behavior is anticipated when 
tij is increased substantially because of the increased re- 
flection at the lead-sample interface). 

IV. TRANSPORT THROUGH STRONGLY 
DISORDERED INTERFACE 

In this section we analyze quantum transport proper- 
ties of a single dirty interface whose dimensionless con- 
ductance G/Gq is much smaller than the number of con- 
ducting channels N c h (i.e., much smaller than the conduc- 
tance of a ballistic conductor of the same cross section). 
For practical purposes, interface can be defined as a any 
scattering region whose thickness is sufficiently shorter 
thanEl Xp (in bulk conductors L 3> A^)- Here we look 
at a geometrical plane of atoms as a model of interface 
in the strict sense. Furthermore, the evolution of trans- 
port from single dirty interface to strongly disordered 
thin slab (i.e., "thin Anderson insulators", since stacking 
together enough finite-size interfaces, with the disorder 
strength chose here, would lead to a bulk Anderson in- 
sulator having exponentially small conductance). Such 
thin slabs are a more likely element in experimental cir- 
cuits!^] These problems are not only conceptual, namely 
to understand the difference between the transport in 
bulk conductors and interfaces, but it has been brought 
about by the understanding of crucial effects the interface 
scattering can have in the CPP transport experimentsE2l 
on GMR magnetic multilayers (in GMR the added com- 
plication is spin-dependent interface resistance^ which 
dominates the resistance and magnetoresistance for not 
too large layer thicknessesE3). 

The importance of interface scattering in many areas 
of metal and semiconductor physics has been realized 
in a plethora of research papers since the seminal work 
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FIG. 9. Conductance of a single dirty interface (N — 1) 
and thin slabs composed of two (N = 2) or three (TV — 3) 
such interfaces (modeled by the Anderson model with diago- 
nal disorder Wi = 30 on a lattice N x 10 x 10). The calculation 
is performed using different values of the hopping parameter 
th in the attached leads, and the results are averaged over 200 
realizations of disorder. 

of Fuchs.@ They are mainly concerned with the trans- 
port parallel to impenetrable rough interface, while re- 
cent have also risen interest in the transport normal to 
the interface (CPP geometry). Because the nature of 
the transport relaxation, time in inhomogeneous systems 
is not well understood^ the first step is to understand 
properties of a single interface before studying them as 
a part of some more complicated circuit like those in 
Sec. 0. For example, the properties of a single interface 
cannot be described in terms of the Boltzmann conduc- 
tivity gb = ne 2 T /m, i.e., using the elastic mean free path 
I = vpT (or transport mean free time r) familiar from 
the bulk metallic conductors whose conductivity is dom- 
inated by the semiclassical effects. 

Lacking enough experimental information, the simple 
theoretical models for the interface effects on electron 
propagation assume diffuse scattering at interdiffuscd 
atoms or interfacial roughnessE^I (these free electron the- 
ories omit the complex electronic band structure of tran- 
sition metals appearing in realistic GMR samples). Even 
disorder-free interface can have a nonzero resistance, be- 
cause of mismatch of the crystal potential and band 
structures .□ Here we are interested only in generic prop- 
erties of mesoscopic transport through interfaces that do 
not depend on material-specific details.cl Therefore, the 
interface roughness is modeled here by the short range 
random scattering potentialu generated by the impuri- 
ties located on the sites of a square lattice 1 x x N z 
and with strong disorder Wj = 30 in the corresponding 
TBH Eq. (jl|). The bulk conductor composed of such 
interfaces (stacked in parallel and coupled with nearest 
neighbor hopping i) is an Anderson insulator, because 
all states are localized already fora W c ~ 16.5. To com- 
pute the CPP transport properties using the mesoscopic 



FIG. 10. Distribution function P(T) of the transmission 
eigenvalues (at half-filling Ef = 0) for a single disordered 
interface N — 1, and thin slabs composed of two TV = 2 
or three N = 3 such interfaces, modeled by the Ander- 
son model with Wi = 30 on a cubic lattice TV x N y x N z . 
The disorder-average is taken over an ensemble of 1000 
conductors. The analytical functions plotted are: (a) 
P(T) = {G/2Gq)1/{T 3/2 VT^ T), an d (b) Dorokhov distri- 
bution P(T) = (G/2Gq) i/(tvT=T). 



Kubo formula ([L4D , the interface is placed in our stan- 
dard computational setup between the two semi-infinite 
disorder-free leads. Thus, the conductance is computed 
for an atomic monolayer of a disordered material inside 
an infinite clean sample of a finite cross section shown 
in Fig. [l]. Also computed is the conductance of a thin 
slab composed of two (i.e., 3D conductor modeled on the 
lattice 2 x N y x N z ) and three sheets (3 x N y x N z ) of 
the same disordered material, as shown in Fig. ||. The 
microscopic origin of different terms contributing to the 
interface resistance was traced back to both specular and 
diffusive scattering in the plane (where diffusive scatter- 
ing can even open additional channels for electron trans- 
port, thereby increasing the conductance) .□ 

Following the discussion in Sec. Ill, the influence of 
the leads on the interface conductance is checked by us- 
ing two different hopping integrals <l (compare to Fig. [?]). 
Here the analysis based on comparison of relevant energy 
scales does not work (i.e., one can not use bulk material 
concepts, like £rh)- Plausibly, I find that leads affect the 
conductance of the interface much more than the conduc- 
tance of a bulk disordered conductor (with similar value 
of disorder- averaged conductance). 

Mesoscopic transport methods give the possibility 
not only to compute the conductance, but also to use 
the picture of conducting channels and their quantum- 
mechanical transmission properties. From a pragmati- 
cal point of view, one does not need these fully quan- 
tum techniques to study the transport in macroscopic 
conductors which are usually dominated by semiclassical 
physics. Nevertheless, the study of the distribution of 
transmission probabilities, which requires phase-coherent 
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FIG. 11. The disorder-averaged (over 200 configurations) 
conductance of a mesoscopic multilayer, composed of dirty 
interfaces and clean bulk conductors, modeled on a lat- 
tice 17 x 10 x 10. The results are obtained from: the 
mesoscopic Kubo formula ([L4]) applied to the whole multi- 
layer (solid line); and using the semiclassical resistor model, 
Grm = (2Ri — Rqpc) 1 (dotted line), where individual resis- 
tances are computed also from Eq. (114) and summed accord- 
ingly (the meaning of Rqpc is explained in Fig. |l5| ). 



transport, enhances our insight into the conduction pro- 
cesses in electronic systems. The diagonalization of a 
Hermitian matrix tt^, which defines the conductance 
through the Landauer formula, 



2e 2 2e 2 ^ 

n—1 



t = 2V-ImS L G\ 



N 



TmEf 



(28a) 



(28b) 



introduces a set of transmission eigenvalues < T n < 1 
for each realization of disorder. The incident flux concen- 
trated in channel \p) will give the wave function in the 
opposite lead ^2 q t pq \q), where t is the transmission ma- 
trix. Here G r 1N is the matrix block, which connects lay- 
ers 1 and N of the sample, of the full retarded real-space 
Green function matrix (15). The distribution function of 
T n is defined as 



P(T) = (J25(T-T n ) 



(29) 



Using P(T), the disorder-average value of any quantity 
that can be written down-in the form of linear statistics 
a(T) can be expressed asE3 



FIG. 12. The disorder-averaged (over 200 configurations) 
conductance of a mesoscopic multilayer, composed of strongly 
disordered interfaces and (quasi) ballistic bulk conductors 
(e.g., £ ~ 9 at half-filling for W — 2), modeled on a lattice 
17 x 10 x 10. The results are obtained from the mesoscopic 
Kubo formula ( |l4| ) applied to the whole multilayer as a single 
quantum-coherent conductor. 



(A) = (jTa(T n )\ = J dTa 



(T)P(T). 



(30) 



For example, the conductance is described by a linear 
statistics g(T) = T. Contrary to the naive expecta- 
tions, stemming from comparison with Drude-Boltzmann 
conductance G ~ GqN^I/ L, that transmission of every 
channel is T ~ £/L in the metallic diffusive (GqN^ 
G ^> Gq) bulk conductor, it was shown by DorokhovE3 
that in a homogeneous multichannel wire geometry 



Pd(T) = 



1 



G 



2G Q TyT^T' 



cosh 



-2 (Nch 

V 9 



< T < 1. 



(31) 



The cutoffl at small T is such that J* dT P(T) = N ch 
ensuring that averages of the first and higher order mo- 
ments of T are not affected (on the proviso ^> 
G/Gq). The Dorokhov distribution is universal — it de- 
pends only on the disorder-averaged conductance G and 
not on the details of disorder, dimension, shape of the 
sample, carrier density, spatial resistivity distribution 
and other sample-specific properties.E3 Thus, P{T) is bi- 
modal distribution function meaning that: most of T n 
are either T n ~ ( "closed" channels) or T n ~ 1 ( "open" 
channels). This has important consequences when calcu- 
lating linear statistics other than the conductance since 
we can get conductance (first moment of the distribu- 
tion) without really knowing the details of P(T) (e.g., 
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FIG. 13. The disorder-averaged (over 200 configurations) 
conductance of a mesoscopic multilayer, composed of strongly 
disordered interfaces and weakly disordered bulk conductors, 
modeled on a lattice 17 x 10 x 10 (£ « 4a at half-filling for 
W = 3). The results are obtained from: (a) mesoscopic Kubo 
formula ( |l4| ) evaluated for the whole multilayer; following the 
resistor model Eq. @, (b) G' KM = (2R A + Rb + 2R I y 1 , and 
(c) Grm = (2R A + Rb + 2Ri - 4Rqpc) _1 . 



FIG. 14. The disorder-averaged (over 200 configurations) 
conductance of a mesoscopic multilayer, composed of strongly 
disordered interfaces and weakly disordered bulk conduc- 
tors, modeled on a lattice 17 x 10 x 10. The results are 
obtained from: (a) mesoscopic Kubo formula (jl4|) evalu- 
ated for the whole multilayer; following the resistor model 
Eq. @, (b) G' RM = (2R A + Rb + 2R1)- 1 , and (c) 
Grm = (2Ra + Rb+ 2Rt - 4Rqpc) _1 - 



the higher moments are probed in the caseL^I of shot-noise 
power, Andreev conductance of metal/supercondcutor 
interface, etc.). 

The simple counting of the number of T n in each bin 
along the interval [0, 1] gives allows us to obtain the dis- 
tribution function P(T) (in this procedure the delta func- 
tion in ( p9| ) is effectively broadened into a box function 
S(x) equal to one inside each bin). Figure plots P(T) 
for the dirty interface and the two thin slabs. The re- 
sult is compared to Pd (T) of Eq. ( ]3l| ) and the one which 
describes the numerical data 



1 



G 



2G Q T 3 / 2 VT^T' 



(32) 



This formula is, up to a numerical factor, rthe same as 
the analytical prediction of Schep and Bauera for a single 
dirty interface 

Thus, interfaces belong to a universality class different 
from that of diffusive bulk conductors characterized by 
the Dorokhov distribution Pd(T). However, it seems that 
distribution is valid even for thin slabs whose thickness 
L is grater than Xp, but is smaller than the localization 
radius (which in our case can be estimated in the pedes- 
trian way as the thickness at which conductance van- 
ishes). This is in compliance with experimental confir- 
mation of Psb(T) in a uniform Nb/AlO^/Nb Josephson 



junction, whose subharmonic gap structure in the I — V 
characteristic is extremely sensitive to the number of con- 
ducting channels at Ep and their transmission probabil- 
itiesEll, since the thickness of the realistic A\O x barrier is 



bigger than Xp in the junction electrodes (Rcf. 20 gives a 
simple derivation of the Schep-Bauer distribution with- 
out using the assumption L <C Xp). In fact, even our 
extremely high disorder input in the standard Anderson 
model generates surface conductance in the band center, 
gj ~ 9 • 10 10 f2 -1 cm~ 2 (for a = 3 A), that is much higher 
than 10 s f2 _1 cm -2 reported for the measured value in 
Ref. [2(] (suggesting that thin slab is more likely to play 
to role of a barrier in these Josephson junction than the 
strict geometrical plane) . While being intriguing concept 
in disorder electron physics, universality can be frustrat- 
ing for the device engineers. Not all features of the trans- 
port through dirty interface are universal 



V. TRANSPORT THROUGH MULTILAYERS 

Armed with the knowledge of transport properties of 
dirty interfaces and metallic homogeneous layers, we can 
now undertake the study of circuits composed of such 
elements. The main feature of our circuits is that they 
are of nanoscale size (~ few Xp) and fully phase-coherent 
(i.e., effectively at zero temperature). The choice of dis- 
order and size of the system is driven by the interest here 
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to explore the patterns of breaking (under the influence 
of quantum effects) of a simple description of the mul- 
tilayer in terms of classical resistor network. Therefore, 
I study such deviations from semiclassics by computing 
the exact zero-temperature conductance of several mul- 
tilayers, as well as of their components, using mesoscopic 
Kubo formula (fl4|). This computation takes into account 
all quantum interference and quantum size effects from 
the onset. In general, resistors can be added according 
to the classical Ohm's law only when their size is larger 
than the dephasing length (in this case the quantum 
features of diffusion can still enter through the trans- 
port properties of individual phase-coherent upijts of size 
Lf/,, example being the weak localization effecto at finite 
temperatures). For example, at high enough tempera- 
tures the system can be partitioned into the cubes of 
a macroscopic size where quantum interference be- 
tween wave functions scattered in different cubes can be 
neglected — this gapkes it possible then to define an in- 
tensive quantity£3 like conductivity, by_aolving a classi- 
cal random resistor network problems.LJ The quantum 
composition law for quantum resistors in a chain is more 
complicated since it contains a phase [Variable depending 
on the characteristics of all scatterers.E30 Thus, because 
of quantum interference effects (like conductance fluctu- 
ationsEa) phase-coherent resistors do not add in series, 
and even for homogeneous samples resistance as a func- 
tion of length is not a self-averaging quantity (therefore 
rcctuiritLg disorder-averaging to restore the Ohmic scal- 
inged'O). In the Landauer-Bttiker scattering formalism, 
resistor network description corresponds to a semiclassi- 
cal concatenation of the scattering matrices of individual 
units. (t in (^8|) is just one block of a full scattering ma- 
trix!^), i.e., the concatenation of "probability scattering 
matrices" of successive disordered regions (obtained by 
replacing each element of the scattering matrix by its 
squared modulusl3). 

The multilayers studied here are composed of three 
bulk conductors joined through two dirty interfaces. The 
whole structure is modeled by the TBH Eq. ( |l6| ) on a 
lattice 17 x 10 x 10 where sixth and twelfth monoatomic 



layers contain the same interface studied in Sec. IV. The 
disorder strength within the plane of the interface is fixed 
to Wi — 30, while disorder inside the bulk layers (com- 
posed of five monolayers) is varied. The disorder strength 
is taken to be the same in two outer layers where dif- 
fusive bulk scattering takes place. This type of multi- 
layer can be viewed as a period of an infinite A/B multi- 
layer:E3 layer of material A on the outside (of resistivity 
Pa and total thickness 2La — 10a, where a is the lat- 
tice spacing) and interlayer material B between the in- 
terfaces (of resistivity ps and thickness Lb = 5a). For 
example, for the chosen disorder strengths W = 3 and 
W = 6 the corresponding bulk material resistivities at 
half-filling are p ~ 130 pilcm and 500 pttcm (assuming 
a = 3A), respectively. I neglect any potential step at 
the interface (caused by the conduction band shift at the 
interfaces). Such multilayers are often described semi- 



14 
=5 12 

io 

8 R 

U 



o 

-a 

o 
U 



80 



60 



40 



n 

o 
p 
a- 
e 

o 

o 

CD 



20 K> 



-2 2 
Fermi Energy 



FIG. 15. Conductance of individual resistors comprising 
the multilayers of Sec. [v|: disordered conductors modeled by 
the Anderson model a lattice 5 x 10 x 10 with diagonal disor- 
der strength, (a) W = 6, and (b) W = 3, and quantum point 
contact conductance (1/-Rqpc) of a clean sample modeled on 
the lattice with the same cross section (i.e., on a lattice sup- 
porting 100 conducting channels, where the maximum of open 
channels with T n = 1 at the Fermi level in such 3D ballistic 
conductorsEa is 68 at \E F \ ~ 0.3). 



i 



classicaly in terms of the resistor modell 



(34) 



where Rrm is the total multilayer resistance, Nb is the 
number of bilayers (I study below just one multilayer 
period Ni, = 1), and Rj is the interface resistance (evalu- 
ated in Sec. IV in Fig. ^|). Thus, the resistor model treats 
both bulk and interface resistances as semiclassical ele- 
ments of a circuit in which resistors add Ohmically in 
a series. Nevertheless, quantum diffusion can be impor- 
tant inside each individual resistor element, as discussed 
above.ES From the measurement of -Rrm as a function 
of the layer thickness, the bulk and interface resistances 
can be extracted experimentally. When quantum inter- 
ference effects become important in the CPP transport, 
this picture breaks down. 

The quantum (i.e., zero-temperature) conductance 
of the multilayers is computed using a standard tool 
throughout the paper — mesoscopic Kubo formula of 
Eq. (|l4|) where whole multilayer is treated as a single 
phase-coherent unit. This approach intrinsically takes 
into account finite-size effects in the problemjij as em- 
phasized in Sec. ||, as well as all single-particle quantum 
interference effects. In all calculations the hopping inte- 
grals throughout the disordered sample and in the leads 
are the same (£l = tc = t). This mean s that no ad- 
ditional scattering, discussed-in Sec. Ill), is generated 
at the lead-sample interfaced The remaining resistance 
can sometimes be interpreted as a series addition of two 
quantum resistors. Namely, the semiclassical limit of the 
two-probe Landauer or Kubo formula for conductance, 
obtained e.g., from the stationary-phase approximatiorO 
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FIG. 16. The disorder-averaged (over 200 configurations) 
conductance of a mesoscopic multilayer, composed of strongly 
disordered interfaces and weakly disordered bulk conductors, 
modeled on a lattice 17 x 10 x 10. The results are ob- 
tained from: (a) mesoscopic Kubo formula (|ll|) evaluated for 
the whole multilayer; and (b) following the resistor model 
Eq. @, G RM = {2R A + Rb- 2R Q p C y 1 - The relative error 
[G — Gbm]/G is plotted as (c). 



of the Green function expression 
of the sample, leads to 



|) for the transmission 



(G>- 



-Rqpc 



'S' 



(35) 



song scattering on lm- 

„2j 



which should be valid for not too i 

purities. The "contact" resistancet-3 -Rqpc = 7r/t/e 2 AT c h is 
nonzero even for ballistic conductor because a finite cross 
section can carry only finite currents (the voltage drop 
in this case occurs at the lead-reservoir interface). Sinii-. 
lar interpretation was given for the interface resistance,c3 
and is expected to be valid also for the interfaces embed- 
ded into a multilayer if there is no coherent scattering 
between adjacent interfaces (e.g., when such effects are 
destroyed by sufficiently strong scattering in the bulk). 
These formulas naturally lead to the resistor model where 
different interfaces and bulk layers contribute to the CPP 
resistance as resistors in series. It is often assumed that 
-Rqpc can be neglected when compared to usually much 
higher resistance of a diffusive sample, thus making the 
choice of the two-probe (or some other) ,-aeometry just 
the matter of computational conveniences 

Since resistor model is expected to be retrieved in the 
limit of completely diffusive .scattering in the bulk and 
lack of phase coherence,HE3 in the quest for substan- 
tial quantum effects induced deviations from this pic- 
ture I start from the opposite limit: a multilayer com- 
posed of clean (ballistic) layers and disorder introduced 
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FIG. 17. The same calculation as in Fig. (the same la- 
bels apply), but for a slightly different multilayer with in- 
creased disorder strength in the two outer layers. 



on the interfaces, W A = W B = and Wj = 30. The 
disorder-averaged result is plotted in Fig. [ll] for two dif- 
ferent thicknesses of the interlayer. Even after disorder- 
averaging, the conductance oscillates as a function of Ep. 
This clearly quantum effect is a consequence of the size 
quantization caused by a coherent interference of elec- 
trons reflected back and fort at the strongly disordered 
interface. After adding the disorder into the outer lay- 
ers, the effect persists albeit with a smaller amplitude 
of oscillations (Fig. |l2|). However, while this is plausi- 
ble because of the interlayer being composed of only few 
atomic monolayers^ (i.e., its length is of the order of \f), 
it is somewhat surprising that oscillations increase with 
increasing the separation between the interfaces. Fur- 
thermore, the oscillating conductance is still quite dif- 
ferent from a pure resonant tunneling conductance peaks 
which would occur at the energies of bound .states if inter- 
faces are replaced by the tunneling barriersc£l (e.g., in our 
model this can be generated by reducing the hopping-Ln- 
tegral i mn between the outer layers and the inter layerEj) . 
Although it is obvious that this phenomenon cannot be 
accounted by the semiclassical resistor model, its answer 
is plotted for the sake of comparison. Also, by adding im- 
purities in the interlayer we can follow the disappearance 
of conductance oscillations (which is akin to the studies of 
disorder effects on conductance quantization in ballistic 
conductorsc3) . The effect has almost vanished at W = 2 
although the mean free path, e.g., around the band cen- 
ter, I « 9 (obtained from Bloch-Boltzmann equation with 
Bom-approximation for the scattering on a single impu- 
ritytil) is still larger than L (i.e., the transport within 
interlayer has not yet reached the limit of fully diffusive 
bulk scattering). 
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To sweep through the transition between fully quan- 
tum and (expected) resistor model semiclassical descrip- 
tion, a disorder is introduced in both layers A and B. 
The results are shown in Figs, [l^ and The disorder 
W = 6 is the strongest one in the Anderson model where 
one can still use the semiclassical picture of transport at 
half-filling (i.e., at E F = we getO I ~ a for W ~ £, 
but localization is postponed to much higher values ofEfl 
W ~ 16.5). The other type of homogenous layer is mod- 
eled with W = 3, which is in the crossover between the 
ballistic and the diffusive transport regime, since £ w 4 at 
i?F = and L^4 = = 5a. The statistical error bars on 
the disorder-averaged conductance (over N con { samples), 
defined as AG = WVax G/N, are smaller than the size of 
the dots (this further clarifies that small conductances of 
the multilayers are not the consequence of strong disor- 
der, but are generated by combining few metallic resis- 
tors and dirty interfaces in a series). The conductances 
of individual layers are plotted in Fig. [l5|. In the first 
case Wa = Wb — 3, the naive application of the re- 
sistor model, where -Rqpc is neglected, leads to G RM 
(Grm = 1/-Rrm is the sum of component resistances) 
being smaller than the quantum G computed for the mul- 
tilayer as a single coherent unit. However, the subtrac- 
tion of four -Rqpc, which brings Grm into the form of 
Eq. (^||), gives conductance higher than Grm > G. Here 
-Rqpc is the resistance of a ballistic conductor ("quan- 
tum point contact") with a cross section A — 100a 2 (see 
Fig. |l5|) . According to Eq. (|3q ) , the resistances of homo- 
geneous layer components of the multilayer, which are 
plotted in Fig. |l5| (or Fig. ^ for the interface alone), are 
the sum of their "intrinsic bulk" resistance^ and -Rqpc- 
Thus, when adding these resistors the final sum should 
contain only one -Rqpc if resistor model is to be com- 
pared with the two-probe resistance of the multilayer. In 
all other cases which include layer with diagonal disor- 
der W = 6, both naive and more careful attempt (where 
(Nr — 1) contact resistances -Rqpc are subtracted from 
the sum of resistance of component layers, where Nr is 
the number of bulk and interface resistances used to ob- 
tain such sum) to use the resistor model give the conduc- 
tances which are greater than the ones computed for the 
multilayer as a single conductor. The explanation which 
invokes only extra scattering on the boundaries between 
different components inside the multilayer, which does 
not work in the first case Wa = Wb = 3, is not the only 
one possible. r— ■ 

It was shown recentlynil that even in weakly disordered 
metals (where Boltzmann resistivity is practically indis- 
tinguishable from the one computed from Kubo formula) 
clear separation of a two-probe resistance into contact 
term and diffusive bulk resistance, as implied by the stan- 
dard arguments of Eq. (|3^), is tempered by localization 
effects [i.e., terms of the order £/L which are neglected 
in Eq. (p5|)1, and therefore possible only in very weakly 
disordered systems. To elucidate the effects responsi- 
ble for the difference between the two ways of evaluating 



the multilayer conductance in Figs. [l3] and |lj, I apply 
the same analysis on a multilayer where dirty interfaces 
are removed (i.e., e m = on the sixth and twelfth plane 
along the x— axis). When Wa Wb such multilayers are 
similar to the junction studied in S ec. I lll l where interface 
is not explicitly modeled as in Sec. |IV|, but appears as a 
boundary between two different homogeneous materials 
brought into a contact. In the opposite case Wa = Wb, 
the sample akin to a homogeneously disordered conduc- 
tor whose conductance was studied in Fig. |t]. The re- 
sults of this calculations, plotted in Figs. [l6] and |l7], hint 
that difference between resistor model conductance Grm 
and conductance G of the multilayer as one indivisible 
quantum-coherent unit is quite minuscule, on the proviso 
that extraneous contact resistance terms are properly 
subtracted. Even though the multilayers are mesoscopic, 
where electrons retaining their phase coherence and are 
subject thereby to quantum interference effects, disorder 
averaging effectively destroys most of random interfer- 
ence terms (a sample self-averages at finite temperatures 
when it can be partitioned into a cubes of size L^, as 
discussed above). 113 The surviving interference terms are 
exemplified by those generated by quantum interference 
effects on special trajectories (i.e., Feynman paths), such, 
as the closed loops responsible for weak localization.Ea 
These nonlocal weak localization effects, arising from the 
interference between the amplitudes along the conjugated 
time-reversed loops which span the whole phase-coherent 
sample, are responsible for G being slightly smaller than 
Grm (localization effects inside the individual resistors, 
summed in the resistor model formula, are taken into ac- 
count even by the semiclassical approach). As expected 
in this picture, the discrepancy is verjz_small at W — 3, 
and is increasing for stronger disorderO W — 6. On the 
other hand, conductance fluctuations sneak in through 
the fluctuations of relative error [G — Grm]/G of the re- 
sistor model as a function of Fermi energy (cf. Figs. [l6| 
and 0). 

VI. CONCLUSIONS 

I have studied different macroscopically inhomoge- 
neous mesoscopic disordered 3D conductors by employing 
both fully quantum description, provided by the meso- 
scopic two-probe Kubo, and semiclassical resistor model. 
Before embarking on the technicalities of computations 
some fundamental issues in the quantum transport the- 
ory have been scrutinized: the relationship between the 
Kubo formula in exact single-particle eigenstate repre- 
sentation and Kubo formula for the finite-size sample in- 
spired by mesoscopic physics; the influence of the macro- 
scopic leads ("measuring apparatus") on the computed 
two-probe conductance; and the transmission and trans- 
port properties of a single dirty interface, which are 
markedly different from a standard notions developed 
for bulk conductors. It was shown, by exact computa- 
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tion on examples of homogenous disordered samples and 
inhomogeneous metal junctions, that mesoscopic Kubo 
formula allows one to obtain the reliable results for the 
static zero-temperature conductance (which includes the 
properties of the attached leads) characterizing nonin- 
teracting quasiparticle transport. On the other hand, 
the evaluation of the traditional Kubo formula, based on 
the exact diagonalization of the respective Hamiltonian, 
ends up with a numerical discrepancy when compared 
to mesoscopic methods applied to homogeneous samples 
(because of the necessity to handle small numerical pa- 
rameters, like the width of Lorentz broadened delta func- 
tion), but qualitative fails in the case of inhomogeneous 
samples. This can be traced back to the conceptual ob- 
stacles in applying methods derived for an infinite sample 
to a finite-size system (even after it is extended through 
periodic boundary conditions) where the problem of dis- 
sipation and effect of the sample boundaries (which is 
crucial for the description of mesoscopic devices) have 
been traditionally bypassed for the sake of computational 
pragmatism. 

The study of a single dirty interface, with specific An- 
derson model type disorder, shows that its transmission 
properties are well-accounted by the Schep-Bauer distri- 
bution, but these formula applies approximately also to 
thin slabs composed of a few of such interfaces (which 
is important for experimental investigations where one 
hardly deals with geometrical planes of theoretical anal- 
ysis). The nanoscale mesoscopic multilayers contain- 
ing such interfaces and ballistic bulk conductors exhibit 
disorder-averaged oscillating conductance as a function of 
Fermi energy. This effect of phase-coherence and quan- 
tum size effect is slowly destroyed upon adding the disor- 
der inside the layers. However, even for diffusive scatter- 
ing in the metallic layer components, smooth disorder- 
averaged multilayer conductances can not be completely 
accounted by the resistor model (which sums layer and in- 
terface resistances as resistors connected in series). Since 
this classical approach works well when interfaces are re- 
moved (up to a tiny localization effects arising from inter- 
ference effects inside the whole sample), we can conclude 
that just a single plane of a strongly disordered material 
is enough to bring new quantum effects into the conduc- 
tance of mesoscopic metallic multilayered structures. 
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